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SECTION  I 


OBJECTIVES  AND  SUMMARY 

The  main  thrust  of  this  work  is  to  predict  the  wake  development  on  theoretical 
grounds,  and  to  test  these  predictions  by  experiments  in  the  Aeronutronic 
Supersonic  Wind  Tunnel  (SWT). 

The  theory  is  in  reality  a  summary  of  the  knowledge  on  linearized  compressible 
wake  behavior  for  both  the  laminar  and  turbulent  cases.  By  considering  a 
fixed  family  of  nozzle  geometries  appropriate  to  GDI  applications,  it  was 
possible  to  derive  general  formulas  for  a  wide  range  of  supply  pressures  and 
temperatures,  exit  Mach  number,  type  of  gas  and  nozzle  wall  temperature.  The 
formulas  give  the  fluid  velocity,  temperature,  etc.,  at  any  point  downstream 
of  the  nozzle  cusp  trailing  edges  in  nondimens ional  form,  and  the  graphs  and 
tables  supplied  allow  one  to  judge  the  effect  of  the  wakes  in  an  arbitrary 
cavity.  A  limitation  inherent  in  this  theory  is  its  weakness  immediately 
downstream  of  the  trailing  edge,  both  due  to  the  nonlinearity  and  the  non¬ 
equilibrium  period  of  flow  adjustment. 

Highly  resolved  experimental  measurements  are  also  being  performed  in  the  SWT 
to  verify  the  theory,  at  least  at  one  set  of  conditions:  Mach  4,  laminar  flow, 
exit  Re  of  100,000,  y  *  1.4  and  variable  Tw  (nozzle  wall  temperature).  A 
representative  set  of  multiple  nozzles  was  installed  in  the  SWT,  capable  of 
internal  cooling  to  simulate  low  Tw/To.  Data  on  the  mean  flow  properties  have 
been  obtained  for  adiabatic  conditions  on  the  wake  centerline;  a  smaller 
volume  of  information  already  exists  on  the  fluctuating  flow  with  and  without 
wall  cooling.  These  data  show  fair  to  ^ood  >,greement  with  the  theory,  and 
they  have  spotlighted  those  areas  of  the  theory  where  obvious  improvements  are 
necessary. 

In  addition  to  these  results,  the  wind  tunnel  diffusion  capability  of  the 
Aeronutronic  setup  was  found  adequate.  Recommendations  for  the  extension  and 
completion  of  the  work  are  indicated. 
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SECTION  II 


EXPERIMENTAL  SETUP 

The  supersonic  portion  of  the  nozzle  contours  is  shaped  according  to  the 
coordinates  of  Table  1.  This  contour  prescription  begins  at  the  throat  and 
terminates  at  the  leading  edge.  At  the  throat  there  is  a  slight  surface  slope 
discontinuity.  The  upstream  end  of  each  cusp  is  a  1  inch  diameter  cylinder 
section  (facing  upstream,  see  figure  1),  and  the  blending  of  this  portion  with 
the  throat  is  accomplished  via  a  0.047  inch  radius  cylinder  section.  Overall 
length  of  each  cusp  is  3  inches,  width  3.10  inches  and  trailing  edge  thickness 
is  0.01  inches. 

Details  of  the  nozzles  appear  on  figures  1  and  2.  Two  half-cusps  were  also 
built,  consisting  of  the  two  symmetric  halves  of  a  regular  cusp,  for  attachment 
on  the  floor  and  ceiling  of  the  test  section.  When  assembled,  the  array  thus 
consists  of  three  complete  two-dimensional  DeLaval  nozzles,  each  of  throat 
height  0.0945  inches,  exit  height  1.113  inches,  width  3.10  inches  and  length 
2.855  inches  (from  throat  to  T.E.)  The  arrangement  is  shovm  on  figure  2  and 
the  nozzle  design  method  is  explained  in  Appendix  A. 

The  array  is  installed  in  the  test  section  of  the  Aeronutronic  Division  Super¬ 
sonic  Wind  Tunnel  (AD/SWT),  which  is  3.369  inches  high  and  3.10  inches  wide 
(figure  3).  The  approach  section  upstream  of  the  nozzles  is  6  inches  long, 
and  is  preceded  by  a  gradual,  shallow  angle  contraction  (about  2  feet  long) 
bringing  the  air  supply  from  the  large  stilling  tank  of  the  facility.  In  this 
fashion  the  supply  air  flow  is  parallel  and  free  of  vortices  and  other  anomalies 
dangerous  to  the  unifonnity  of  the  stream. 

The  rectangular  cross  section  continues  for  about  12  inches  downstream  of  the 
array  T.E.  (this  is  the  "cavity"  proper)  at  which  point  the  air  enters  the 
supersonic  diffuser.  No  effort  was  made  to  change  the  diffuser  from  its  normal, 
fixed  configuration,  which  consists  of  a  two-dimensional  constriction  caused  by 
floor  and  ceiling  "ramps".  However,  a  third  ramp  was  later  added  on  the  side- 
wall  to  constrict  the  throat  further  and  improve  the  unchoked  flow.  Further 
details  on  the  diffuser  type  and  operation  appear  in  paragraph  4.2. 

Initially  optical  glass  sidewalls  were  used,  stretching  from  a  point  6  inches 
upstream  of  the  nozzles  to  the  diffuser  entrance.  This  allowed  an  extensive 
Schlieren  study  of  the  flow  which  proved  very  useful  in  detecting  streamlines 
in  the  approach  section,  evaluating  the  nozzle  flow  and  resolving  early  problems 
of  tunnel  choking.  After  completion  of  the  optical  studies,  lucite  windows 
were  installed,  especially  to  avoid  glass  breakage  occurring  when  the  nozzles 
were  cooled.  The  lucite  windows  were  geometrically  similar  to  the  earlier  glass 
windows,  allowing  the  same  ease  of  visual  inspection  of  the  flow. 

The  method  of  cooling  the  nozzles  to  reduce  the  ratio  Tw/To  and  study  its 
effect,  consists  of  circulating  liquid  nitrogen  (LN2)  through  them.  The  cooling 
channel  is  a  single  hole  running  through  the  nozzle,  coaxial  with  Che  main  struts 
supporting  Che  nozzles  on  the  sidewall.  Cooling  was  confined  only  on  the  two 


TABLE  1 


TABLE  OF  PHYSICAL  COORDINATES  OF  CUSP  NOZZLES 


^  (in.) 

X  (in.) 

p  (in.) 

X 

0 

.0398 

.5142 

.4896 

.1808 

.4089 

.052 

.482 

.1916 

.4034 

.0575 

.4784 

.2029 

.3978 

.0618 

.4759 

.2153 

.3918 

.0646 

.4741 

.2272 

.3862 

.0674 

.4726 

.2398 

.3803 

.0699 

.4711 

.2530 

.3744 

.0734 

.469 

.2658 

.3686 

.0836 

.4627 

.2945 

.3563 

.0873 

.4603 

.3251 

.3437 

.0912 

.4579 

.3582 

.3307 

.0943 

.4562 

.394 

.3172 

.0971 

.4545 

.4326 

.3034 

.100 

.4527 

.4748 

.2889 

.1025 

.4511 

.5215 

.2742 

.1074 

.449 

.5715 

.2592 

.1125 

.4463 

.6274 

.2433 

.117 

.4436 

.6866 

.2274 

.1214 

.4407 

.7525 

.2111 

.1268 

.4376 

.904 

.1768 

.137 

.432 

1.0888 

.141 

.146 

.427 

1.3161 

.1045 

.1529 

.4235 

1.5902 

.0704 

.1616 

.4189 

1.9317 

.0407 

.1704 

.4143 

2.3479 

.0176 

1  'X. 

2.6157 

2.8545 

.0076 

.005 

full  nozzle  cusps,  since  Che  half  nozzles  attached  to  the  tunnel  floor  and 
ceiling  could  not  be  effectively  cooled  because  of  conduction.  Two  thermo¬ 
couples  were  installed  on  one  of  the  cusps  to  monitor  its  temperature,  one  on 
Che  forward  (subsonic)  portion  of  the  cusp  and  the  second  after  the  nozzle 
throat.  The  coolant  was  drawn  from  LN2  "bottles",  whose  capacity  is  adequate 
to  produce  uninterrupted  wind  tunnel  runs  lasting  many  hours. 
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SECTION  III 


INSTRUMENTATION 


3.1  DIAGNOSTICS  OF  TUNNEL  OPERATION  AND  PERFORMANCE 

A  group  of  instruments  is  "fixed"  on  the  wind  tunnel  and  monitors  overall' 
facility  operation,  and  nozzle  and  diffuser  operation  during  the  run.  These 
measure : 

3.1.1  Stagnation  (total  or  "supply")  pressure  Po 

3.1.2  Stagnation  (total  or  "supply")  temperature  To 

3.1.3  Moisture  content  of  flowing  air 

3.1.4  Temperatures  at  other  points  of  the  circuit  (e.g.  ambient, 
diffuser,  etc.) 

3.1.5  Pressure  midway  between  the  supersonic  and  subsonic  diffusers 

3.1.6  Pressure  downstream  of  the  subsonic  diffuser 

3.1.7  Barometric  pressure 

3.2  FLOW-FIELD  INSTRUMENTATION 

A  minimum  of  three  independent  measurements  are  sufficient  to  define  the  flow 
field  at  each  point.  Commonly,  these  consist  of  the  pitot  tube,  the  static- 
pressure  probe  and  the  To  probe.  The  pitot  tubes  used  in  this  work  are  con¬ 
strained  mainly  by  their  size  (e.g.  the  outside  diameter  (O.D.)  of  the  tip) 
which  must  meet  resolution  requirements,  and  their  "hydraulic  volume"  which 
must  meet  time-response  requirements.  Adequate  resolution  can  be  obtained 
if  the  OD  is  less  than  0.1  of  thn  width  of  the  shear  flow  in  question.  For 
almost  the  entire  flow  field  iti  these  experiments,  this  requirement  v;as  sat¬ 
isfied  by  the  0.02  inches  OD  pitot  tubes  used.  Using  thin  tubing  to  connect 
the  pitot  tube  to  its  transducer,  the  time  constant  could  be  held  to  1  second 
or  less.* 

A  variety  of  static  pressure  probes  are  available,  a  typical  one  consisting 
of  a  0.02  inch  tube  tapered  to  a  sharp  point  at  the  upstream  end,  with  four 
0.006  inches  diameter  holes  drilled  circumferencial ly  at  a  distance  of  about 
0.25  inches  from  the  tip.  Viscous  interaction  is  computed  to  be  negligible 
for  this  geometry. 

The  pitot  and  static  pressure  probes  are  connected  pneumatically  to  a  "pressure 
station"  positioned  outside  the  test  section.  This  station  mounts  two  absolute 

*In  continuous  surveys  (i.e.  probe  moving  continuously)  the  adequacy  of  Che 
response  was  judged  by  repeating  the  same  traverse  at  progressively  diminishing 
probe  speeds.  Maximum  permissible  speed  was  judged  to  be  the  one  below  which 
no  "lag"  was  seen  in  the  analog  traces. 


pressure  transducers  (0-2  psia  and  0-5  psia  range),  a  bank  of  Wallace-Tiernan 
absolute  pressure  gages  (0-20,  0-100  and  0-400  mmHg  absolute),  a  digital  volt¬ 
meter  and  associated  conditioning  electronics,  a  network  of  valves  and  manifolds 
and  a  vacuum  pump.  These  self-contained  systems  give  great  flexibility  to  the 
operator  during  the  run.  For  example,  the  pitot  or  static  pressure  can  be 
obtained  directly  as  a  voltage  or  read  on  the  gages,  or  the  transducer  cali¬ 
bration  can  be  checked  quickly.  Usually  the  transducer  in  use  is  calibrated 
just  before  and  after  a  data-taking  run.  The  accuracy  of  the  Wallace-Tiernan 
gages  is  checked  periodically  by  the  Metrology  Laboratory  of  the  Aeronutronic 
Division.  Where  digital  data  is  desired,  the  system  output  is  channeled 
directly  to  the  Supersonic  Wind  Tunnel  Data  Acquisition  System.  When  connected 
to  an  analog  (XY)  recorder,  the  system  can  rapidly  produce  "pitot"  maps  which 
can  reveal  small  details  of  the  flow. 

The  total  temperature  is  measured  locally  with  the  so-called  Tq  probe  which 
is  usually  a  bare  wire  thermocouple.  As  with  all  probes  measuring  recovery 
temperature,  this  probe  is  sensitive  to  the  Mach  and  Reynolds  number;  a 
"recovery  factor  calibration"  is  then  an  important  prerequisite  to  using  this 
probe.  This  sensor  also  has  the  simplest  conditioning  electronics,  requiring 
only  that  one  thermo junction  is  immersed  in  ice. 

Hot  wire  and  film  anemometers  are  used  for  the  detection  and  measurement  of 
turbulence.  No  such  quantitative  data  were  obtained  in  the  present  period, 
and,  thus,  it  is  unnecessary  to  describe  the  complex  hardware  and  software 
needed  for  such  quantitative  measurements.  Qualitative  observations  in  these 
tests  were  made  with  hot  film  anemometers  newly  developed  at  Aeronutronic 
Division.  Useful  bandwidth  for  these  extended  to  beyond  300  KHz, 
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SECTION  IV 


JOURNAL  OF  OPERATING  CHARACTERISTICS  OF  THE  NOZZLE  ARRAY 

4.1  OVERALL  CHARACTER  OF  NOZZLE  FLOW 

Apart  from  minor  corrections  made  in  the  adjustment  of  the  nozzles,  the  array 
generally  operated  well  from  the  start.  Figure  4  shows  a  "pitot  probe  map"  of 
the  flow  obtained  at  z  «  0,  in  which  the  wakes  and  wake  shocks  are  clearly 
seen.  Figure  5  is  the  corresponding  picture  of  the  flow  obtained  by  a  one 
microsecond  Schlieren  mosaic  of  the  flow  in  the  cavity.  In  the  approach  section 
the  only  extraordinary  feature  is  the  apparent  separation  of  the  tunnel  floor 
and  ceiling  boundary  layers  at  the  upstream  ends  of  the  half  cusps.  This  was 
fully  expected  and  does  not  seem  to  influence  the  flow  in  the  upper  and  lower 
nozzles.  In  the  nozzles  themselves  and  beyond,  one  can  see  clearly  the  nozzle 
boundary  layers  and  wakes,  the  boundary  layers  growing  on  the  floor  and  ceiling 
of  the  cavity  (test  section),  the  wake  shocks  emanating  at  the  T.E.  and  the 
diffuser  inlet  shocks. 

4.2  TUNNEL  AND  DIFFUSER  OPERATION 

It  should  be  mentioned  at  the  outset  that  it  was  not  the  objective  of  this 
work  to  optimize  the  diffuser  performance  of  the  tunnel,  but  rather  to  operate 
the  latter  in  the  "unchoked"  mode  so  that  the  growth  of  the  wakes  would  not 
be  disturbed. 

For  the  type  of  flow  pictured  in  figure  5,  conditions  were  Pq  *  754  mmHg  and 
Tg  ■  146°F.  These  conditions  prevailed  for  almost  the  entire  set  of  tests 
done  so  far.  However,  even  when  these  conditions  were  set,  occasional  choking 
of  the  tunnel  was  initially  observed,  pictured  in  figures  6  and  7.  The  former 
shows  a  time  record  of  the  output  of  the  pressure  transducer  sensing  downstream 
of  the  supersonic  diffuser.  When  this  output  is  of  order  2  and  3  millivolts 
the  tunnel  is  unchoked  and  choked,  respectively.  In  this  figure  it  is  shown 
that  the  tunnel  was  unchoked  for  about  1.6  hours  after  start  and  then  lapsed 
into  a  half  hour  period  of  choked  flow.  In  the  choking  cycle  the  shocks  are 
observed  to  move  first  to  a  position  upstream  of  the  supersonic  diffuser  in¬ 
let  and  progress  gradually  to  the  nozzle  exits  (figure  7).  In  the  "unchoking" 
portion  of  the  cycle  the  process  reverses  itself  until  the  shock  system  is 
again  expelled  to  a  position  downstream  of  the  diffuser  throat. 

Since  no  tunnel  control  changes  lay  behind  the  appearance  of  choking,  its 
occurrence  is  inexplainable .  Some  attempts  to  improve  the  recovery  were  made, 
consisting  mainly  by  decreasing  the  diffuser  throat  opening  from  8.3  square 
inches  to  7.3  square  inches.  The  latter  was  the  opening  employed  throughout 
the  tests;  further  efforts  to  improve  recovery  were  abandoned,  since  several 
hours  of  continuous  unchoked  operation  is  now  routinely  easy  to  obtain  in  the 
facility. 

The  diffuser  throat  area  Aq  is  thus  0.71  of  the  test  section  area  A.  The  total 
pressure  drop  across  the  nozzle  diffuser  system,  for  Pq  ■  754  mmHg,  is  694  mmHg, 
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giving  p£  «  60  tmiHg  aC  Che  pump  Inlec.  Thus,  Che  operaclng  pressure  racio  is 
12.5,  as  opposed  co  Che  normal  shock  raCio  of  7.4  (for  Che  design  Mach  4). 

This  corresponds  Co  a  pressure  recovery  of  60%  of  Che  normal  shock  recovery 
which  is  noc  bad  considering  chac  no  efforc  was  made  Co  modify  subscancially 
Che  diffuser  of  Che  facilicy.  Furchermore,  che  nacural  nonuniformicy  of  che 
flow  in  a  "GDL  cavicy"  raises  che  possibiliCy,  as  is  now  widely  underscood, 
chac  che  recovery  in  diffusers  wich  nonuniform  inleC  flows  is  noc  a  simple 
problem. 

4.3  NOZZLE  OPERATION  AND  ALIGNMENT 

The  nozzle  array  operaced  well  from  che  mechanical  viei^oinc.  Some  crouble 
was  aC  firsc  encouncered  wich  che  mechod  of  accaching  che  upper  and  lower  half¬ 
nozzles  (figure  3)  on  Che  floor  and  ceiling  of  che  channel.  Uhcil  chese  were 
finally  epoxied  scrongly  in  place,  Cheir  0.005  inch  chick  T.E.  cended  Co 
decach  from  che  surface  and  fluccer  wich  a  "whiscle"  clearly  audible  Co  che 
cunnel  operacor. 

The  alignmenc  of  Che  nozzles  was  found  Co  be  cricical  in  order  Co  obtain  an 
absolucely  uniform  flow  ac  cheir  exic.  The  procedure  followed  was  Co  inscall 
che  nozzles  and  Chen  survey  downscream  of  chem  wich  che  picoc  cube.  Asymmecries 
found  by  chese  piCoc  craverses  led  Co  a  readjuscmenc  of  che  nozzles,  and  che 
cycle  was  repeaced  uncil  uniform  flow  was  obcained.  Criceria  for  "unifonnicy" 
were  che  symmecry  of  che  wake  profiles  and  of  che  wake  shocks  and  cheir  locacion, 
che  idencicy  of  che  growch  races  of  che  cwo  wakes,  laceral  (spanwise)  uniformity, 
ecc.  In  che  beginning,  when  each  sidewall  was  a  single  piece  of  opcical  glass, 
Chis  iceracive  process  was  very  cumbersome  because  probe  changes  also  caused 
nozzle  alignmenc  loss  (since  che  nozzles  were  discurbed  when  che  cesc  sec  cion 
probe  was  accessed.)  Window  breakage  during  che  nozzle  cooling  tescs  prompced 
replacemenc  of  Che  glass  sidewall  on  each  side  wich  a  lucice  Cwo  piece  window, 
che  ups  cream  half  of  which  supporced  che  nozzles.  This  was  che  "final"  nozzle 
supporc  arrangemenc  which  allows  once  and  for  all  fixing  of  che  nozzles  in  place 
(during  each  Cesc  encry),  visibilicy  of  che  encire  flow  field*  and  easy  access 
Co  che  Cesc  secCion  probe. 

4.4  COOLING  SYSTEM  OPERATION 

The  nozzles  were  cooled  by  circulating  through  the  liquid  nitrogen  (LN2) 
supplied  under  pressure  from  bottles  (each  with  160  liter  capacity.)  Two 
thermocouples,  one  ac  Che  stagnation  point  of  Che  cusp  and  another  located 
0.5  inches  downscream  of  the  throat  and  flush  wich  che  surface,  indicated  che 
nozzle  temperature  Tw  (where  Tw  is  usually  the  average  of  the  cwo  readings.)  Low 
temperature  would  be  achieved  very  rapidly,  one  minute  being  typically  che  time 
it  cook  Co  lower  che  Tw  Co  che  desired  value.  Intermediate  Tw  could  be  obcained 
at  will,  and  held  steady  by  throttling  Che  bottle  valve.  Typical  operating 
conditions  were: 


*Usually,  no  Schlieren  observations  are  made  with  chis  arrangemenc  because  of 
-che  poor  opcical  properties  of  lucice. 
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Stagnation  temperature  To  »  146°F  (606°R) 
Minimum  Tw  achieved  -  -100°F  (360°R) 
Adiabatic  nozzle  temperature  ■  90°F  (550°R) 
Adiabatic  Tw/To  ■  0.91 
Minimum  Tw/To  ■  0.59 


Note  that  the  minimum  Tw/To,  although  not  excessive,  is  ample  for  the  purpose 
of  checking  the  Tw/To  effect  on  wake  growth  which,  by  reference  1,  is  kno^m 
to  be  important.  For  example,  for  the  conditions  of  this  experiment,  (M*4, 

Y  **  1.4)  reference  1  predicts  that  at  a  distance  of  10  nozzle  exits  the  temp¬ 
erature  defect  on  the  wake  centerline  is  1.4  for  the  laminar  adiabatic  case 
(Tw/To  ■  0.9),  and  1.0  for  Tw/To  =  0.6.  Further  decreases  in  the  attainable 
Tw/To  are  thought  possible  by  minor  rework  of  the  cooling  ducts  inside  the 
nozzle  cusps.  (Paragraph  5.3  for  the  effect  of  Tw  on  transition). 

During  the  cooling  tests  the  flowing  air  dev^joint  was  kept,  as  always,  at  or 
below  the  -35°F  mark  by  periodic  reactivation  of  the  tunnel's  silica  gel  air 
dryer.  Even  so,  a  fine  layer  of  frost  appeared  on  the  upstream  side  of  the 
cooled  cusps,  but  not  on  the  nozzle  surfaces.  This  is  an  Indication  that  both 
friction  and  heat  transfer  is  more  pronounced  on  the  latter  surfaces,  as 
assumed  in  reference  1.  No  adverse  effect  of  this  frost  layer  was  noted  on 
the  flow. 


9 


SECTION  V 


EXPERIMENTAL  RESULTS 


5 . 1  NOZZLE  FLOW 

A  profile  of  Mach  number  along  2  ■  0,  x  -  0  (che  nozale  exit  plane)  is  shown 
on  figure  S.  It  is  seen  that  the  design  exit  Mach  of  4.0  has  been  well 
achieved,  although  on  the  nozzle  axes  (y  "  0,  y  ■  i  h^ )  there  is  a  slight 
increase  to  M  ■  4.10.  A  boundary  layer  profile  at  the  exit  is  shown  on 
figure  9,  as  obtained  with  the  pitot  cube.  It  is  seen  that  the  edge  5  of  che 
boundary  layer  at  this  position  (T.E.)  is  almost  exactly  equal  to  0.10  inches. 
A  representative  theoretical  prediction,  utilizing  che  formula  of  Low  (ref¬ 
erence  2)  for  laminar  flow  5  on  a  flat  place  with  pressure  gradient. 
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gives  6  ■  0.075  inches.  This  does  not  compare  badly  with  the  experiment, 
considering  che  theoretical  approximations  involved  and  che  fact  chat  6  as 
deduced  from  figure  9  will  be  larger  than  the  actual  6  due  to  Prandtl  number 
effects . 

This  comparison  indicates  strongly  chat  che  nozzle  boundary  layers  are  laminar. 
A  second  indication  is  the  shape  of  che  pitot  profile  of  figure  9  which  is 
completely  in  character  with  laminar  flow.  A  third  indication,  equally  as 
strong,  is  chat  the  value  of  Che  unit  Reynolds  number  at  che  exit  is  100,000 
inch's  so  chat  che  wetted- length  Reynolds  number  Re^  ■  280,000.  The  latter 
is  an  order  of  magnitude  below  the  transition  Reynolds  number  for  flat  plate 
boundary  layers,  as  explained  in  Appendix  B.  Using  the  results  of  this 
appendix,  it  is  estimated  chat  pq  should  be  increased  to  about  100  psia  in 
order  to  achieve  natural  transition  on  che  nozzle  surfaces  in  the  present 
experimenta’  setup. 

Figure  10  shows  lateral  static  pressure  surveys  taken  for  the  middle  nozzle  at 
its  exit  plane.  Here  the  sensor  was  fixed  first  at  y  ■  0  (che  vertical  center 
of  che  nozzle  exit  plane  at  x  ■  0)  and  then  traversed  from  z  ■  0.6  inches  to 
-0.6  inches,  which  includes  a  little  less  than  half  che  span  (figure  2).  By 
thus  traversing  che  sensor  parallel  to  che  two-dimensional  nozzle  throat, 
possible  departures  from  flow  two-dimensionality  can  be  revealed.  The  absolute 
level  of  che  static  pressure  for  this  po  (*754  mmHg  abs)  at  M  ■  4  is  also  clear 
in  che  figure.  Following  che  traverse  for  che  middle  nozzle  (labeled  "y»0"  in 
the  middle  of  che  figure),  che  traverses  for  the  lower  (at  left,  for  y  *  -0.97") 
and  the  upper  (y  ■  +0.97")  nozzle  are  shown.  A  small  disturbance  at  the  edge 
of  che  y  ■  -0.97  inches  trace  is  che  only  exception  to  an  otherwise  very  smooth 
two-dimensional  flow  picture.  This  disturbance  is  due  to  che  "throat  wave" 
emanating  from  che  intersection  of  the  middle  nozzle  throat  and  che  wind  tunnel 
sidewall.  It  is  not  clear  from  the  work  done  so  far  chat  such  throat  waves 
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emanate  from  all  such  intersections,  but  the  disturbance  is  in  anyway  small, 
as  the  figure  indicates. 

Figure  11  shows  the  measured  Mach  number  along  the  centerline  of  each  of  the 
three  nozzles  as  far  downstream  as  6hi  .  There  is  a  small  systematic  decrease 
of  M,  but  more  serious  is  the  appearance  of  disturbances  which  cause  local  M 
variations  of  order  10%.  At  least  part  of  this  problem  is  suspected  to  be 
due  to  waves  emanating  from  the  wind  tunnel  sidewalls  (floor  and  celling). 
Figure  12  shows  additional  measurements  of  Mach  number  in  and  downstream  of 
the  nozzles.  The  upper  part  of  the  figure  shows  data  along  the  three  center- 
lines  (y  ■  z  ■  0)  and  the  bottom  part  shows  profiles  along  y  both  along 
the  vertical  line  x  >  z  >  0  and  at  2.4hi  downstream  of  the  nozzle  exit. 

In  summary,  the  comments  made  in  this  Section,  and  an  Inspection  of  the  ac¬ 
companying  figures,  reveal  the  following  for  the  nozzle  flow: 

1.  The  design  Mach  number  goal  was  achieved  within  a  few  percent. 

2.  Good  two-dimensionality  and  repeatability  is  achieved. 

3.  Laminar  boundary  layers  exist  on  the  nozzles. 

4.  Within  certain  small  deviations,  supersonic  flow  is  maintained  in 
the  cavity  and  at  the  nozzle  exit  Mach  number  for  a  distance  of  at 
least  6  times  the  nozzle  exit  heights. 

5.2  QUANTITATIVE  WAKE  MEASUREMENTS 

I 

In  the  present  test  series,  quantitative  wake  measurements  were  done  within 
the  following  matrix: 

Pq  ■  754  mmHg  abs 
To  -  146°F 
y  “iO.Sblhi 
z  »  0 

0  <  X  <  6hj 

Re'  ■  97,200  inch"^  (nominal) 

CO 

Mo,  ■  4  (nominal) 

Tw/To  -  0.92 

which  indicates  that  the  centerline  ("axis"  or  "centerplane")  of  the  upper 
and  lower  wakes  were  probed  at  near  adiabatic  (uncooled)  conditions. 
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Figure  13  shows  Che  measured  variation  of  static  pressure  p  and  local  total 
temperature  To.  The  former  generally  increases  away  from  the  T.E.,  an  indica¬ 
tion  of  Che  expected  flow  "slowdown”  as  the  diffuser  entrance  is  approached. 
"Bumps"  in  this  distribution  are  not  bad,  considering  that  p  is  an  extremely 
sensitive  parameter.  For  example,  it  has  been  calculated  that  Che  pressure 
increase  shown  in  1.8  <  x/hi  <2.7  for  the  lower  wake  would  be  generated  by  a 
"shock  generator"  of  wedge  angle  ~  1.5°.  Such  small  flow  deflections  are 
typical  of  what  is  expected  for  Che  present  setup;  e.g.,  at  the  backward 
facing  step  of  the  nozzle  T.E.  It  is  doubtful  that  any  substantial  improve¬ 
ment  in  the  pressure  distribution  could  be  made  even  if  more  effort  were  to 
be  invested  into  nozzle  design  improvement.  The  wake  To  variation,  shown  on 
the  same  figure,  shows  the  characteristic  drop  below  the  To„  level  and  the 
gradual  Increase  doxms cream. 

The  important  question  in  this  work  is  whether  the  wake  behavior  agrees  with 
the  theoretical  predictions,  now  collected  in  reference  1.  This  reference 
indicates  that  Che  velocity,  temperature,  density,  and  total  temperature  . 
defects  w,  t,  r  and  9  respectively,  defined  by 


w 


u„  -  u(0) 


C  T(0) . -  T, 

T_ 


p«  "  p(0) 
r  3 - - 

0  00 


To(0)  -  To„ 


should  be  unique  functions  of  the  coordinate  Ax/hi  ,  where  Ax  s  x  -  x^  is  Che 
distance  from  the  virtual  wake  origin  Xg  (  and  h^  is  the  nozzle  exit  height). 
Thus,  to  compare  Che  values  of  w,  r,  t,  etc.,  measured  in  this  experiment 
with  the  theoretical  predictions  it  is  important,  first,  to  see  where  the 
origin  Xq  lies  according  to  the  experiment. 

5.2.1  DETERMINATION  OF  THE  VIRTUAL  ORIGIN  FROM  THE  DATA 

The  procedure  to  be  followed  is:  (a)  a  test  is  first  made  to  see  if  the 
relation  . 

w,  c  (x  -  (1) 
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holds,  (b)  points  lying  close  Co  Che  T.E.  and  noc  obeying  che  above  relation 
are  excluded  from  partaking  in  che  finding  of  the  origin  Xq,  (c)  a  computer 
least- squares  fit  of  Che  remainder  of  che  points  is  made  and  Xq  is  found. 

Figure  14  shows  Che  w,  r,  t  plots  versus  Ax/hj  .  For  Che  moment  we  will  use 
this  figure  to  test  che  validity  of  che  theoretical  prediction  (1).  Indeed, 
this  relation  is  seen  to  be  achieved  if  che  first  six  points  (at  the  lowest 
Ax/hi  )  are  excluded.  Leaving  figure  14  aside  for  a  moment,  make  a  second 
test  of  the  relation  (1)  by  plotting,  more  conveniently,  the  inverse  square 
of  w  and  t,  for  both  wakes,  on  figure  15,  versus  x/hj  .  Open  symbols  in  this 
plot  are  che  "excluded"  six  points  in  each  case.  All  other  points  (solid 
symbols)  were  computer  least-squares  fitted  into  a  straight  line.  The  inter¬ 
sections  of  these  lines  with  the  abscissa  should  give  che  virtual  origin  x^. 
There  are  differences  for  che  upper  and  lower  wakes  as  follows: 

Upper  wake: 

Velocity-defect  determination:  Xq  ■  -3.6hj 

Temperature-defect  determination:  Xq  ■  -3.8hj 

Average:  -3.7hj 

Lower  wake: 

Velocity-defect  determination:  Xg  ■  -0.6hi 
Temperature-defect  determination:  x^  »  -1.2hj 

Average:  -0.9hi 

Average  for  both  wakes:  Xg  ■  -2.3hj 

The  large  spread  obtained  between  che  two  wakes  is  less  than  satisfactory, 
although  it  is  clear  that  xo  lies  upstresn  of  che  T.E.  As  a  temporary  measure, 
therefore,  proceed  by  adopting  che  theoretically  predicted  value  Xg  ■  2.5h 
which  is  close  to  che  "grand  average." 

5.2.2  VARIATION  OF  THE  MEASURED  DEFECTS 

The  data  on  the  measured  values  of  w,  t,  r  and  S  are  shown  on  figures  14  and 
16.  These  plots  show  information  additional  to  chat  of  figure  15  in  that  che 
r  and  ^  are  now  included,  che  coordinates  are  logarithmic,  and  use  has  been 
made  of  che  fact  chat  che  effective  origin  lies  at  Xg  ■  2.5hi .  In  this  way, 
for  example,  che  first  point  taken,  at  x  ■  0.061hi  ,  is  shown  on  figure  14  to 
lie  at  Ax/hi  ■  2.5  +  0,061  ■  2.561.  In  this  way,  che  curves  for  the  theoretical 
predictions  from  reference  1  have  been  drawn  in  for  comparison,  since  in  this 
theory  it  was  also  assumed  that  Xg  ■  -2.5hi  for  the  conditions  M,,,  y>  etc., 
of  the  present  experiment.  Furthermore,  two  secs  of  curves  are  shown  for  the 
theory  (solid  and  dashed);  che  dashed  curves  correspond  to  predictions  for  the 
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design  conditions  M®  =»  4,0,  y  “  Tw  =  To.  The  solid  curves  show  the 

predictions  for  the  flow  conditions  actually  obtained  (which  were  slightly 
different  from  design)  as  follows:  Moo  *  3.75*,  Tw  ~  0.9  To. 

The  main  results  from  figures  14  and  16  are  twofold:  first,  there  is  a  clear 
agreement  with  the  one-half  power  wake  decay  law,  and  the  location  and  extent 
of  the  nonequilibrium  region  is  also  clear;  theta  findings  were  already 
gleaned  from  figure  15.  The  second  important  result  is  that  there  is  some 
disagreement  with  the  theoretical  predictions.  Specifically,  it  appears  that 
the  measured  defects  are  about  50%  lower  than  the  predictions,  on  the  average. 
From  this  finding  it  is  tentatively  concluded  that  the  theory  is  not  especially 
accurate.  An  examination  of  the  basis  of  the  theory  (reference  1)  shows  that 
the  latter  includes  linearization  and  a  few  other  subtle  assumptions  which  are 
normally  "safe".  The  present  results  indicate  the  need  for  improving  the 
theory . 

It  should  be  noted  that  the  observed  disagreement  is  in  a  very  sensitive  type 
of  experimental  quantity,  the  "defect"  (see  the  earlier  definition  of  w,  t, 
r  and  9).  In  numerical  (dimensional)  terms,  the  fluid  axis  velocity,  temper¬ 
ature,  etc.,  as  measured,  agrees  with  the  predictions  closely,  (See  discussion 
to  this  effect  on  p.  33  of  reference  1). 

An  effort  was  made  to  find  additional  experimental  data  to  compare  with  the 
theory  of  reference  1.  A  great  number  of  GDL  setups  have  already  been  studied, 
but  unfortunately  most  studies  are  limited  to  Schlieren  photography.  Of  those 
which  are  not,  the  data  of  Director  (reference  3)  provide  numerical  information 
on  wake  decay  similar  to  those  obtained  here  (Director's  data  concern  only  the 
density  defect  detected  by  interferometry).  These  results  are  shown  on  figure 
17  compared  with  the  theory;  the  agreement  is  good,  and  is  significant  because 
his  wake  was  turbulent. 

5.3  FLUCTUATION  MEASUREMENTS  WITH  AND  WITHOUT  COOLING 

Quantitative  measurements  of  the  "turbulence"  of  the  wakes  were  made  only  in 
a  preliminary  fashion,  partly  because  the  precise  measurement  of  turbulence 
is  too  big  a  task  to  occur  so  early  in  this  research,  and  partly  because  the 
wake  was  already  suspected  to  be  laminar.  A  thin  film  anemometer  was  chosen, 
chiefly  because  its  sturdiness  makes  it  more  suitable  for  exploratory  work. 
Traverses  of  the  rms  signal  were  made  versus  y  for  a  number  of  constant  x 
and  constant  z. 

Figure  18  shows  typical  results  when  the  nozzle  cusps  were  uncooled  (e.g.. 

To  =  146°F,  Po  »  754  mmHg,  Tw  =  95°F).  At  midspan  (z=0)  and  a  representative 
axial  distance  x  *  hj  the  figure  shows  the  wideband  sensor  output,  character¬ 
ized  of  the  two  double  peaked  wake  profiles.  The  signal  between  the  two 

*As  already  explained,  the  nozzle  exit  Mach  number  Mo,  was  found  a  little  dif¬ 
ferent,  4.15  vs.  4.0,  at  the  center  of  the  nozzle.  On  the  wake  axis,  however, 
the  Mach  number  started  out  at  about  4.0  and  decreased  to  about  3.75  down¬ 
stream.  Therefore,  it  was  felt  that  a  good  alternate  estimate  of  the  pre¬ 
dictions  would  be  to  use  M  =  3.75. 
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wakes  is  almost  exactly  equal  to  the  electronic  noise,  i.e.,  there  is  no 
sensible  signal  in  the  space  bounded  by  the  two  wakes.  The  wake  signals  are 
higher  than  this  noise  signal  by  a  factor  of  about  4.  The  question  posed 
here  is  this:  if  this  wake  signal  results  from  wake  turbulence,  then  the 
claim  of  a  laminar  wake  for  these  conditions  is  contradicted.  The  dilemna 
is  resolved  by  studying  the  sensor  output  spectra  shown  on  figure  18.  It  is 
seen  that  the  energy-carrying  spectrum  lies  below  25  KHZ,  On  the  other  hand, 
the  "wake  frequency"  computed  for  these  wakes  on  the  basis  of  their  width 
(~  0.2"  =  0.5  cm)  and  the  edge  velocity  (  ~  700  m/sec)  is  140  KHZ,  so  that 
the  spectrum  should  extend  beyond  1  MHZ  if  the  wake  was  turbulent.*  Thus,  the 
wake  is  indeed  laminar,  and  the  unanswered  question  remaining  is  the  origin  of 
the  observed  fluctuations  in  this  laminar  wake.  A  study  of  earlier  experiments 
(reference  4)  is  that  fluctuations  in  the  laminar  wake  exist  long  before  tran¬ 
sition  sets  in. 

Similar  measurements  were  made  with  the  nozzle  cusps  cooled  to  T„  ~  -80°F 
(corresponding  to  T^/Tq  ~  0.63).  Here  it  was  found  that  the  wake  thickened 
by  a  factor  of  about  2  from  x*0  to  x*6hi  where  the  fluctuation  spectrum  cas¬ 
caded  to  about  300  KHZ.  It,  therefore,  appears  that  the  cooling  action  "tripped" 
the  wake  into  turbulence.  This  confirms  the  well  known  phenomenon  of  wake  de- 
stabilization  by  cooling  (reference  5).  Future  measurements  are  badly  needed  to 
pinpoint  the  accelerated  transition  quantitatively. 


*A  practical  rule  for  estimating  the  upper  "active"  frequency  f  of  a 
turbulent  shear  flow  is  f  =  lOue/b,  where  Ug  is  the  edge  velocity  and 
b  is  the  flow  width. 
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SECTION  VI 


INTERIM  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  WORK 

On  the  basis  of  the  work  done  to  date,  the  following  conclusions  can  be  drawn 
about  the  specific  experiment  performed  in  the  AD/ SWT: 


(1)  The  wind  tunnel  was  successfully  converted  into  a  facility  for  testing 
GDL  nozzle  arrays. 

(2)  A  nozzle  array  was  installed  and  tested,  fulfilling  its  design  requirements 
at  Mach  4.0  and  continuous  flow. 

(3)  Although  no  design  provisions  were  made  to  "tailor"  the  nozzle-wall  inter¬ 
section,  good  quality  two-dimensional  flow  was  achieved  to  x  *  6h^  . 

(4)  Running  time  is  arbitrary  (e.g.,  several  hours)  with  unchoked  flow  without 
drastic  changes  in  the  diffuser  configuration. 

(5)  The  cooling  system  can  be  further  improved,  although  lonj  'uration  flows 
at  T^  ^  0.6  Tq  can  be  already  achieved  at  present. 

In  addition,  certain  general  conclusions  of  importance  to  all  zzle  cusp  wake 
research  can  be  drawn: 

(6)  The  quantitative  data  taken  so  far  are  in  fair-to-good  agree  unt  with  a 
rudimentary  theory  of  wake  flows  (reference  1). 

(7)  The  virtual  wake  origin  for  nozzle  cusp  wakes  appears  to  be  located  at  the 
nozzle  throat. 

(8)  Nozzle  wall  cooling  moves  turbulent  tr='nsition  upstream  in  the  wake. 

(9)  The  trailing  edge  nonequilibrium  region  occupies  a  length  equal  to  about 
two  nozzle  exit  heights,  beyond  which  wake  similarity  is  achieved. 

In  order  of  importance,  the  following  casks  are  recommended  for  the  immediate 
future : 

(10)  Measurements  of  the  mean  (average)  properties  must  be  continued  and  com¬ 
pleted.  Specifically,  the  following  must  be  done: 

10.1  Data  must  be  taken  all  the  way  to  the  diffuser  entrance  (x  ~  12h,  ) . 

10.2  Radial  (lateral)  profile  data  must  be  obtained. 

10.3  A  better  determination  of  the  virtual  origin  is  possible,  from  the  data 


10.4  A  complete  set  of  data  must  be  obtained  for  the  cold  wall  case,  and 
compared  with  the  theory  of  reference  1. 

(11)  Turbulent  wakes  must  be  obtained  and  measured.  Transition  by  artificial 
roughening  of  the  nozzle  surface  should  be  possible  to  obtain. 

(12)  The  cause  of  the  observed  fluctuations  in  the  laminar  wake  (mentioned  in 
paragraph  5.3)  must  be  found  by  further  hot  wire  anemometer  measurements. 

(13)  For  the  turbulent  wake,  hot  wire  anemometer  measurements  must  be  combined 
with  a  set  of  turbulence-magnitude  predictions,  especially  for  the  cold  wall 
case,  to  see  if  it  is  possible  to  predict  the  fluctuations  for  nozzle  cusp 
wakes  (in  the  same  fashion  as  the  mean  flow  properties  have  been  predicted  by 
reference  1). 
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Figure  5.  Spark-Schlieren  photos  (shovm  at  two  different  cont 
of  the  flow.  Note  boundary  layer  separation  in  the 
chamber  (at  right)  and  diffuser  inlet  shocks  at  lef 
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Figure  13.  Stacie  pressure  and  total  temperature  as  measured  along  the 
axes  of  the  two  wakes  (z  -  o)- 
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Figure  14.  The  temperature  (t),  density  (r)  and  velocity  (w)  defects 
along  the  two  wakes  compared  with  theoretical  predictions- 
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Figure  17.  The  data  of  Director  compared  with  the  Aeronutronic  theory. 


film  of  the  r  01  a  fluctuatlona  In  the  flow  (at  left) 
onpanied  by  a  apectruai  aurvey  at  right. 


APPENDIX  A 


DESIGN  OF  THE  GDL  NOZZLES 


A.l  DESIGN  SOURCE 

The  nozzles  were  designed  according  Co  Che  mechod  of  Cwo-dimensional  characcer- 
isCics  as  puc  forch  in  Shames  and  Seashore  (reference  6).  This  reference 
designs  for  Che  shorcesc,  highesc  expansion  nozzles  possible.  Originally,  che 
"scraighc forward"  mechod  of  reference  6,  p.  6-8,  was  accempced,  buc  Che  wall 
concour  resulced  in  a  final  area  abouc  107.  higher  chan  chac  given  for  isencropic 
expansion.  Therefore,  che  mechod  finally  followed  is  che  "averaging"  mechod. 

A. 2  AVERAGING  METHOD 

Figure  A.l  shows  che  averaging  mechod.  The  line  ab  is  drawn  from  the  inicial 
wall  poinc  a,  ac  an  angle  PoinC  b  is  che  incerseccion  of  chis  line  with 
che  characteristic  leaving  Che  kernel  at  poinc  K  ac  an  angle  \.  From  che 
midpoint  c  of  line  ab,  a  line  is  drawn  ac  angle  ty' ,  which  is  che  wall  inclina¬ 
tion  following  inclination  in  che  reference  6  cable.  This  intercepts  che 
previous  characteristic  from  K  at  B,  and  che  new  characteristic  from  K'  (at 
angle  \')  ac  d.  The  new  wall  point  is  B,  and  che  new  '\nidpoint"  is  e.  The 
coordinates  of  the  kernel  points  K  and  che  inicial  wall  point  a  are  given  in 
Che  tables  of  reference  6  in  terms  of  che  "half  throat  height"  A|;/2,  as  are 
Che  angles  nt  and  \.  In  che  algorithm,  a  and  B  are  called  "wall  points",  b 
and  d  "far  points"  and  c  and  e  "midpoints". 

A. 3  CALCULATION  ALGORITHM 

A  Cable  was  sec  up  as  follows: 


TABLE  A- I 

FUNCTIONAL  DIAGRAM  FOR  COMPUTING  NOZZLE  COORDINATES 


\ 

0  -  - 


\  a 


Can\ij 


cana  Bj^  B^ 


Y  X 

1  0  wall  poinc 

-  -  far  point 

-  -  midpoint 


0.01  - 
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using  Che  reference  6  cables  (for  Mf  ■  4,  see  p.  13-14)  Co  lisc  Y"  chrough  can  ^ 
(che  qualicies  and  lisced  here  are  equal  Co  X/Af/2  and  Y/A^/Z  of  chose 
Cables).  The  lasc  chree  columns  were  lefc  blank.  The  coordinaces  of  Che  firsc 
wall  poinC,  1  and  0  were  also  puc  in.  Referring  co  figure  A.l  and  Table  A-I 
now,  Che  equaCion  for  Che  lines  ab  and  Kb  are: 


y  ■  X  CancTa  +  (Xa  “  ^a  ^  ®wa  (A.l) 

y  ■  X  tanXjj  +  (yj^  -  Xj^  canXj^)  3  x  CanX^  +  (A.2) 

Since  Che  inCercepC  of  lines  passing  chrough  poincs  K  depend  only  on  Che 
Xr,  yjj,  Xfj  aC  ChaC  poinc,  che  column  B^^  in  Table  A-I  was  also  filled  in 
advance : 

\  %  canXk  (A.3) 


The  sCeps  were  as  follows: 

(1)  The  inCercepc  of  line  ab  was  deCermined  from  ics  def  ini  cion: 


Bw  s  Xa  -  ^a  '=^°“a 

(■  1  -  0  ■  1  in  chis  case  only)  (A, 4) 

(solid  lines  in  Table  A-I) 


(2)  Point  b  was  determined  from  che  incersection  on  lines  ab,  kb,  i.e,,  by 
solving  simulcaneously  eqs.  (A.l)  and  (A.2): 


tarofa  -  CanXfj 

yb  - - 

CanXj^  -  tanc^a 


Bw  -  Bk 

Xb  - - 

CanXj^  -  canofg 

(docced  lines  in  Table  A-I) 

(3)  The  midpoinc  c  was  decerrained  from: 

Xa  +  Xb 

Xc  -  -  ,  y^,  -  - 

2  2 


(A.5) 


(A.6) 


(A.7) 
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(in  this  first  case. 


0  +• 


X 


(double  arrows  on  Table  A- 1)  ^  Une 

K  toire  X.  y/  and  the  tan  a  from  the  next  lin 
L)  A  new  is  determined  by  using  Jc. 
of  Table  A-I  and  apply  equation  (a.4). 

.  yc  -  x^  tan  a 

(dashed  lines  in  Table  A-I) 

,,  ac  B  IS  tiov  Obtained  by  applying  fo^las  <a.5) 

(1  Values  enclosed  by  wavy  lines  in  Table  A-I. 

«;tb  tbe  new  wall  point  tbns  bnown.  tbe  ptocednte  again  begins  with 

Step  !• 

A. 4  RESULTS 

,.4.i  NONDIMENSIONAL  COORDINATES  WITHOUT  VISCOUS  CORRECTION 

This  procedure  yielded  the  “fr  tL^tilt^dJc^arfd  ToTisel'' 

algorithn  gave  an  inafarea  ratio  Af/A^  -  10.80  instead 

tropic  expansion  to  Mach  .  ,  discrepancy  was  corrected  by  displacing 

t  tfe^'no^tu  l^t  S:ril:o  Jp^sed  because  we  felt  it  a  bettet  design 
pcactice;  this  will  be  discussed  futther. 

d  ■55(A./2)  =  27.5At  -  2.565Af,  which  is  within 

5rt^:ft^rrtb:ore"faf  "Sd0tt‘n5tale'’  length  of  I.BObAg. 
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TABLE  A- II 


INVISCID  NOZZLE  COORDINATES  (non-dimensional) 


X 

y 

At/2 

At/2 

0 

1 

.768 

1.496 

1.002 

1.647 

1.109 

1.715 

1.191 

1.767 

1.246 

1.800 

1.300 

1.834 

1.347 

1.864 

1.416 

1.908 

1.611 

2.037 

1.684 

2.087 

1.759 

2.135 

1.817 

2.171 

1.872 

2.206 

1.928 

2.244 

1.975 

2.276 

2.071 

2.316 

2.168 

2.375 

2.256 

2.429 

2.34 

2.489 

2.445 

2.552 

2.641 

2.692 

2.815 

2.793 

2.947 

2.842 

3.115 

2.937 

3.284 

3.031 

X 

y 

At/2 

3.484 

3.141 

3.693 

3.253 

3.910 

3.367 

4.15 

3.49 

4.378 

3.604 

4.621 

3.723 

4.875 

3.845 

5.123 

3.962 

5.676 

4.213 

6.265 

4.468 

6.902 

4.731 

7.593 

5.002 

8.336 

5.279 

9.15 

5.565 

10.05 

5.865 

11.012 

6.166 

12.09 

6.482 

13.23 

6.797 

14.5 

7.118 

17.42 

7.788 

20.98 

8.478 

25.36 

9.169 

30.64 

9.79* 

37.22 

10.295* 

45.24 

10.635* 

50.4 

10.71** 

55 

10.72** 

*CorrecCed  for  over expans ion 

**Added  to  ensure  proper  isentropic  expansion 
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X.  the  distance  from  the  throat.  Is  used  to  compute  Re,,  but  Re'  -  , 

value  comes  from  local  conditions. 

Other  conditions  for  the  calculation  were  T„  -  -  100°F,  -  730  mmHg. 

The  resulting  physical  thriisplaceme^It  Sickness  at  the 

text).  To  compute  them  fectionnhe^ght  H  -  3.3688  in.  Between 

exit,  6^  -  0-0^95  inches,  h,  -  .  1.1229". 

centerlines  of  the  cusp  trailing  edges,  tH  3 

Allouln,  0.01"  UP  thlchness.  the  physical  half-helRht  alloued  is 


h  1.1229-0.01 


0.5565" 


h  -  1.113" 

The.,  if  the  throat  height  is  A.  end  the  Elou  height  at  the  exit  is  h. 

h  -  hf  +  2S*  ■  W.^2ic  +  26*  ■  b 

or 

10.72A^  +  .099  »  1.113 


A  -  0.0946" 

Knowing  A^,  the  coordinates  x,  y  Jg^^errthen  found  from  6*.  Finally, 

?Sn'o:idfnrrf-  tht  ‘c:sf  crterpLne  to  the  cusp  surface  (Table  1) 
was  found  as  an  aid  to  the  machinist. 
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'ure  A.L  Schematic  of  nozzle  design  method  and  coordinate 


APPENDIX  B 


TRANSITION  LOCATION  IN  THE  GDL  NOZZLE  FLOW 
B.l.  PURPOSE  AND  METHOD 

The  objeccive  was  Co  predict  Che  location  of  trtasition  on  and  downstream  of 
the  GDL  nozzles  as  presently  under  test.  Such  an  estimate  can  be  made 
easily  if  the  nozzle  boundary  layer  is  approximated  to  a  flat  place  boundary 
layer;  i.e.  if  each  point  on  Che  nozzle  is  assumed  to  lie  on  a  flat  plate  of 
edge  M  and  Re'  as  found  on  that  point.  Standard  flat  plate  transition  correla¬ 
tions  from  experiments  supplied  another  necessary  datum,  as  is  explained  below. 

B.2.  TRANSITION  CORRELATION 

A  flat  plate  correlation  of  laminar  boundary  layer  transition,  for  adiabatic 
walls,  was  obtained  from  Coles  (reference  7)  in  the  usual  form  (,<»,  Rex)  in 
Che  plane  of  Mach  M  and  wetted  Reynolds  Re^. 

B.3.  NOZZLE  FLOW  TRAJECTORY 

A  pair  of  Mg(*Ma,)  and  Re'  (unit  Reynolds  number)  obtains  at  each  point  on  Che 
nozzle,  as  well  as  a  pair  Mg  and  Re^^  s  xRe',  where  x  is  the  distance  from  the 
throat  parallel  to  the  synmetry  axis.  The  Mg,  Re^^  curve  is  Che  "nozzle  tra¬ 
jectory"  and  transition  is  assumed  to  occur  where  this  curVe  intercepts  the 
transition  correlation. 

The  procedure  for  plotting  the  Mg,  Rex  curve  was  computerized  co  apply 
more  generally.  The  first  input  was  the  known  variation  of  nozzle  surface 
contour: 


y  ■  y(x),  y*  s  y(x-o)  ■  y(M»i)  (B.l) 

where  y*  (the  throat  height)  is  another  known  number.  The  contour  (B.l)  was 
curve-fitted  as  follows: 

y  -  E  Cn  n-0,1,2 _  (B.2) 

n 

where  C  are  coefficients  found  for  each  nozzle  under  study,  and  y/y*  or  A/ A* 
was  formed.  To  avoid  using  the  implicit  isentropic  relation  for  the  area 
ratio)  where  scars  refer  to  the  conditions  at  the  throat) 


A  A 
A**  A* 


(M) 


Che  relation  was  inverted  and  curve-fitted  into 


(B.3) 


^  2  = 
n  ^ 


(isentropic  area  ratio) 


(E.4) 
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At  this  point,  it  was  decided  to  further  account  for  the  ratio  of  specific 
heats  Yj  so  that  the  final  result  for  local  M  on  the  nozzle  as  a  function 
of  area  ratio  A/ A*  and  y  was 

f  -  0.7111Y  +  0.9956  ^2  2) 

“  "  1.9956  -  0.7111Y 

where  f  is  the  (iaverted  )isentropic  M  -  A/ A*  relation  from  (B.4). 

To  proceed,  the  stagnation  viscosity 

(Tg  +  460)^^^ 

u,  =  1.09  X  10“5  T  +  658'. 6 

was  computed,  where  T  is  the  stagnation  temperature  in  °F.  Finally,  the 
unit  Reynolds  number  for  stagnation  pressure  Pq  at  each  M  was  found  from 


where  the  isentropic  temperature  ratio  T/T^  is  computed  at  M  and  where 
is  now  in  degrees  Rankine. 

B.4  FIRST  TEST  CASE:  TRANSITION  IN  THE  MACH  3.0  NOZZLE  OF  THE  SWT 

Using  equation  (B.7),  the  trajectory  Rex,  M  was  computed  for  the  Mach  3  nozzle. 
Since  M  ■  M(x),  the  Rex,  x  trajectory  was  also  computed  and  is  shown  on  figure 
B.l,  superposed  on  Cole's  transition  criterion.  Transition  is  predicted  when 
Pq  ■  730  mmHg  absolute,  but  not  when  Pg  ■  200  mmHg,  a  verified  fact. 

Figure  B.2  is  an  extension  of  figure  B.l,  plotting  Pg,  x  instead  of  Rex,  x. 

The  transition  prediction  is  generally  good  at  large  x,  judging  from  Schlieren 
data,  but  poorer  at  small  x.  Near  the  nozzle  throat,  especially,  transition 
is  observed  to  occur  much  earlier  than  predicted. 

B.5  SECOND  TEST  CASE:  MACH  4  GDL  NOZZLES 

Figure  B.3  shows  the  nozzle  trajectory  as  computed  from  equation  (B.7)  for  the 
highest  available  P^  for  continuous  flow.  The  highest  Re^^  attained  at  the 
nozzle  exit  is  jurtc  short  of  3  x  lO^.  The  small  Rex  is  mainly  due  to  the  short 
nozzle  length  and,  to  a  lesser  extent,  due  to  the  decreased  unit  Re'  for  the 
Mach  4  operation. 
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Figure  B.3  predicts  chat  the  GDL  nozzle  surface  boundary  layers  will  be 
laminar,  assuming  smooth  adiabatic  wall  and  the  assumptions  leading  to  Che 
computation  of  Re', 

B.6  TRANSITION  IN  THE  WAKES 

Very  slender  bodies  are  known  Co  produce  unstable  wakes  at  supersonic  speeds, 
with  transition  Reynolds  number  occurring  around  60,000.  Since  this  Reynolds 
number  is  based  on  flow  properties  external  to  the  wake  and  Che  distance  from 
Che  body,  transition  in  the  GDL  nozzle  wakes  is  expected  to  occur  within 
1  cm  of  the  cusp  T.E.  Generally,  however,  this  subject  is  quite  complex  and 
will  be  approached  again  in  a  later  publication. 
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Prediction  of  transition  in  the  nozzle  laminar 
boundary  layer  (SWT  Mach  3  nozzle). 


I 


I 


X  (INCH) 

REYNOLDS  NO.  VArTaTION  -  GDL  MACH  4  NOZZLES 


Figure  B.3.  Prediction  of  transition  in  the  nozzle  boundary 
layer  for  the  Mach  4  GDL  nozzles.  In  this  case 
transition  does  not  obtain. 
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